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We revisit Fn{Q^) and Fp^{Q'^), P = 'K,r),r)' , making use of the local-duality (LD) version of 

QCD sum rules. We give arguments, that the LD sum rule provides reliable predictions for these 

form factors at Q^ > 5 — 6 GeV^, the accuracy of the method increasing with Q^ in this region. For 

the pion elastic form factor, the well-measured data at small Q give a hint that the LD limit may 

be reached already at relatively low values of momentum transfers, Q^ « 4 — 8 GeV^; we therefore 

^vj ' conclude that large deviations from LD in the region Q^ = 20 — 50 GeV^ seem very unlikely. The 

data on the {ri,Tj') — > 77* form factors meet the expectations from the LD model. However, the 

f~*) ■ BaBar results for the vr" — >■ 77* form factor imply a violation of LD growing with Q^ even at 

^SJ ' Q^ ~ 40 GeV^, at odds with the ri, 77' case and with the general properties expected for the LD sum 
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INTRODUCTION 



(— I ' In spite of the long history of theoretical investigations of the pion, its properties are still not fully understood. 

O ' At asymptotically large values of the momentum transfer, Q^ — ^ cx), the pion elastic and the 7r7 transition form 

I ■ factors obey perturbative QCD (pQCD) factorization theorems [l|, 01 

^ : F.{Q^) -> 8^a,(g2)/2/g2, F,^(Q2) ^ ^/2^/Q^ /„ = 130 MeV. (1.1) 

\ Subleading logarithmic and power corrections modify the behaviour (jl.ip at large but finite Q^ . 

In early applications of QCD to the pion elastic form factor, Ft^{Q^), one hoped that power corrections vanish 
rather fast with Q^; however, later investigations revealed that nonperturbative power corrections dominate the form 
factor Ft^{Q'^) up to relatively high Q^ « 10-20 GeV^. This picture has arisen from different approaches [2-la|- I^ ^^^ 
\j^ found that even at Q^ as large as Q^ = 20 GeV^ the 0{\) term provides about half of the form factor; the pQCD 
fS| formula based on factorization starts to work well only at Q^ > 50-100 GeV^. 

^J ■ For the pion-photon transition form factor, Ft^jIQ"^), axial anomaly [3, [l3| fixes its value at Q^ = 0. Brodsky 

^\ \ and Lepage proposed a simple interpolating formula between the known value of the form factor at Q^ = and its 
^vq ' asymptotic behaviour (11.11) 
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2 ' which was believed to describe well the transition form factor in a broad range of Q^. 

■ - - Surprisingly, some of the recent studies of the pion elastic form factor in the region Q^ ~ 5 — 50 GeV^ [lll - [l4l | 

reported much larger valus of the pion elastic form factor than expected before (see Fig. [T]); the BaBar result on 
the TT -> 77* transition form factor [16[ imply a strong violation of pQCD factorization in the region of Q^ up to 40 



GeV2. 

In this paper, we revisit F^{Q'^) and Fp^{Q^) making use of the local-duality (LD) version [13, [3 of QCD sum 
rules [l9| . Our main emphasis is on the analysis of the expected accuracy of this powerful (although approximate) 
method for the calculation of hadron form factors. 

A local-duality sum rule is a dispersive three-point sum rule at r = (i.e., infinitely large Borel mass parameter). 
In this case all power corrections vanish and all details of the non-perturbative dynamics are hidden in one quantity 
- the Q^-dependent effective continuum threshold. Implementing duality in the standard way, i.e., as the low-energy 
cut in the dispersion representations for the form factors, the sum rules relate the pion form factors to the low-energy 
region of Feynman diagrams of pQCD: 



F^ iQ)^j^ J dsi I ds2 Apcrt(si,S2,Q'), F^}/ [Q') ^ — I dsap,,tis,Q'). (1.3) 
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Fig. 1: Some predictions for the pion elastic form factor -/'^(Q^) — lower solid (black) line: BLM'2008 [g|, upper solid (green) 
line: BT'2008 [H, dashed (magenta) hne: GR'2008 [li|, and dash-dotted (blue) line: BPS'2Q09 [H, Eq. (4.11b)] — vs. experi- 
ment UM- 



Here Aport(si, S2, Q^) is the double spectral density of the (AVA) 3-point function; whereas crport(s, Q^) is the single 
spectral density of the (AVV) 3-point function. These quantities are calculated as power series in as'- 

Apert(si,S2,Q') = A|,clt(si,S2,g') + a,A|,^i,t(si,S2,Q2)-|-0(a2), ap„.t(s, Q^) = 4ort(s, Q') + ^(a^). (1.4) 
The one-loop spectral densities Ap^j.^ and a^J^^ are well-known |2(]| - |22| . The two-loop contribution A^J^^ has been 



calculated in [23|; the two- loop 0{as) correction to iTpert was found to be zero [2J]. Higher-order radiative corrections 
are unknown. 

If one knows the effective thresholds Soff (Q^) and Scff (Q^), Eqs. (|1.3p provide the form factors. However, finding a 
reliable criterion for fixing the thresholds is a very subtle and difficult problem investigated in great detail in [23 • 

Due to specific properties of the spectral functions at large values of the momentum transfer, the LD form factors 
(|1.3p obey the factorization theorems (II. ip as soon as the effective thresholds satisfy the following relations: 



Soff 



{Q' 



00) = Scff((9^ ^ 00) = 47r2/^- 



(1.5) 



For finite Q^, however, the effective thresholds Scs{Q^) and Soff(Q^) depend on Q^ and differ from each other [2f 
The "conventional LD model" arises if one assumes (jl.Sp for all "not too small" values of Q^ [ITj] : 



scff(g') - seff(g2) = An^fl 



(1.6) 



Obviously, the LD model (II. 6p for the effective continuum thresholds is an approximation which does not take into 
account details of the confinement dynamics. The only property of theory relevant for this model is factorization of 
hard form factors at large values of the monetum transfers. Because of the approximate character of the predictions 
from LD sum rules, it is important to understand the expected accuracy of the form factors obtained by this method. 
Quantum-mechanical potential models provide a possibility to probe this accuracy: one calculates the exact form 
factors by making use of the solutions of the Schrodinger equation and confronts these results with the application of 
LD sum rules in quantum mechanics j27H29l |. 

This paper is organized as follows: In the next section, we recall details of LD sum rules in QCD and of the LD 
model for form factors. Section 3 studies the accuracy of the LD model for elastic and transition form factors in 
a quantum-mechanical potential model. The pion elastic form factor is discussed in Section 4 and the P — >■ 77* 
transition from factors are considered in Section 5. Section 6 gives our conclusions. Appendix A provides technical 
details of perturbative two-loop calculations in non-relativistic field theory. 



2. LOCAL-DUALITY MODEL FOR FORM FACTORS IN QCD 



The sum-rule calculations of the hadron form factors are based on the OPE for the relevant correlation functions, 
which contain the contributions of ground-state hadrons of interest. In order to extract the form factors of pseudoscalar 
mesons, one analyses the {AVA) and (AVV) correlators, A being the axial and V the vector currents. 



A. The three-point function (AVA) and the pion elastic form factor 

The basic objects for the extraction of the pion decay constant and the elastic form factor are the two- and three- 
point correlation functions 



n(/) - J d'xe^P^{0\T[jix)jHO)]\0), 
r{pl,plq') = fd^x^d^X2e'^P^^^-P^^^^Q\T[j{x^)J{Q)jHx2)]\0), q^Pi~P2, Q''^-q^; (2.1) 



where |0) is the physical QCD vacuum, which differs from pQCD vacuum; properties of the physical QCD vacuum are 
characterized by the condensates [19| . j is shorthand for the interpolating axial current jsq of the positively charged 
pion, 

(0|j5a(0)K(p))=ipaA. (2.2) 



J labels the electromagnetic current J^, and for brevity we omit all Lorentz indices. In QCD, the correlators (|2.ip 
can be found by applying OPE. Instead of discussing the Green functions (12. ip in Minkowski space, it is convenient to 
study the time-evolution operators in Euclidean space, which arise upon performing the Borel transformation p" — > r 
to a parameter r related to Euclidean time. The Borel image of the two-point correlator n(p^) is 

oo 
noPE(T) = / ds e^"^ /Opert(s) + ncond(T), /9pcrt(s) = po{s) + Os pi{s) + O(a^), (2.3) 



with spectral densities Pi{s) related to perturbative two-point graphs, and nonperturbative power corrections Ilcond (■'')• 
At hadron level, insertion of intermediate hadron states casts the Borel-transformed two-point correlator into the form 

n(r) = /J e"™' "• + excited states. (2.4) 

In this expression for n(T), the first term on the right-hand side constitutes the pion contribution. Applying the double 
Borel transform p^ 2 ~^ ■''/2 to the three-point correlator r{pl,p2,q^) results, at QCD level, in 

00 00 

roPE('r,Q^) == / /dsids2expf ^^-^ r j Apcrt(si, S2, Q^) + rcond(T, Q^), 



Aport(si, S2, Q^) = Ao(si, S2, g^) + Us Ai(si, §2, Q^) + O(a^), (2.5) 

where Aport(si, S2, Q^) is the double spectral density of the three-point graphs of perturbation theory and rcond(''', Q^) 
labels the power corrections. Inserting intermediate hadron states yields, for the hadron-level expression for T{t, Q^), 

r(T, Q2) ^ F^iQ^) fl e""'^ + excited states. (2.6) 

Quark-hadron duality assumes that above effective continuum thresholds Scu the excited-state contributions are dual 
to the high-energy regions of the perturbative graphs. In this case, the relevant sum rules read in the chiral limit 

Sctt(-r) / 2\ - 2 

L = ds e Ppcrt(s) + ^^^ T + — T -I- • • • , (2.7) 





FAQ'')fl = J J dsids2Ap„.t(si,S2,Q')exp 



Si + S2 



^{^,^^_I^,.^rS + Q^r)+-. (2.8) 

As a consequence of the use of local condensates, the right-hand side of (12. 8p involves polynomials in Q^ and therefore 
increases with Q^, whereas the form factor F.^{Q'^) on the left-hand side should decrease with Q^. Hence, at large Q^ 



the sum rule (|2.8I) . with its truncated series of power corrections, cannot be directly used. There are essentially two 
ways for considering the region of large Q^. 

One remedy is the resummation of power corrections: the resummed power corrections decrease with increasing Q^. 
This may be achieved by the introduction of nonlocal condensates [3l| in a, however, model-dependent manner 32] . 

Another option is to fix the Borcl parameter t to the value r = 0, thus arriving at a local-duality sum rule 17'. 18i|. 
Therein all power corrections vanish and the remaining perturbative term decreases with Q^. In the LD limit, one 
finds 

f. = J d.ppe.t(.) = g (l + ^) + 0{al) , (2.9) 





FAQ') I; =11 dsids2Apert(si,S2,Q'). (2.10) 



The spectral densities /Opert(s) and Apert(si, S2, Q'^) are calculable by perturbation theory. Hence, by fixing Soff and 
Scff(Q^), it is straightforward to extract the pion's decay constant /tt and form factor Ft^{Q'^). 

Noteworthy, the effective and the physical thresholds are different quantities: The latter is a constant determined by 
the masses of the hadron states. The effective thresholds ScR and Scg are parameters of the sum-rule method related to 
the specific realization of quark-hadron duality; in general, they are not constant but depend on external kinematical 
variables (25ll26J. 

Let us recall the important properties of the spectral densities on the right-hand sides of (|2.3|) and (|2.5p : For Q^ — >• 0, 
the Ward identity relates the spectral densities pi{s) and Ai(si, S2, Q^) of two- and three-point functions to each other: 

lim A,(si,S2,Q^)= P*(si) (5(51-52), z = 0,l,.... (2.11) 

For Q^ — > oo and si^2 kept fixed, explicit calculations [S] yield 

1 Q 

lim Ao(si,S2,(5^) ex -— , lim Ai(si, S2, Q^) = 7^ Po(si) po(s2). (2.12) 

For the pion form factor F^(Q^) on the left-hand side of (|2.8p two exact properties are known, namely, its normalization 
condition related to current conservation, requiring -FV(O) = 1, and the factorization theorem ()1.1|) . Obviously, for 

Seff (Q' ^ 0) = /; , Seff (g' ^ C^) = SLD = ^1^^ fl (2.13) 

1 -I- as(0)/7r 

the form factor F^((5^) extracted from the LD sum rule (|2.10p satisfies both of these rigorous constraints. At small Q^, 
we assume a freezing of as{Q^) at the level 0.3, as is frequently done. At the intermediate Q^ , the effective threshold 
depends on Q"^ and the pion form factor obtained from the LD sum rule depends on the details of Seff(Q^)- The 
"conventional LD model" assumes that reasonable predictions for the elastic form factor at not too small values of 
Q^ may be obtained by setting Seff(Q^) = 47r^/^, see Eq. (|1.6p . 

Although one has to invoke assumptions on behaviour of the effective threshold at the intermediate Q^, some 
features of the pion form factor turn out to be largely independent of this assumption. So, let us look more carefully 
what is in fact conjectured in the LD sum rule and what may be predicted by this approach. 

The sum rule (|2.10|) for the pion form factor relies on two ingredients: first, on the rigorous calculation of the 
spectral densities of the perturbative-QCD diagrams (recall that power corrections vanish in the LD limit r = 0); 
second, on the assumption of quark-hadron duality, which claims that the contributions of the hadronic continuum 
states may be well described by the diagrams of perturbation theory above some effective threshold Soff. Thus, the 
only — although really essential — unknown ingredient of the LD sum rule for the pion elastic from factor is this 
effective continuum threshold Scff(Q^). Let us emphasize that, since the 0(1) and 0{as) contributions to the pion form 
factor are governed by one and the same effective threshold Scff (Q^), the relative weights of these contributions may 
be predicted. Their ratio fI^ '{Q'^)/Ft^ (Q^) turns out to be relatively stable with respect to ScS and may therefore 
be calculated relatively accurately (see Fig. [5]). 

Quark-hadron duality implies that the effective threshold (although being a function of Q^) always — i.e., also for 
large Q^ — stays in a region close to 1 GeV^. Moreover, in order to satisfy the QCD factorization theorem for the 
0{as) contribution to the form factor, the effective threshold should behave like Soff(Q^) —> 47r^/^ for Q^ — >■ oo. This 
requirement has immediate consequences for the large-Q^ behaviour of different contributions to the pion form factor: 
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Fig. 2: Ratio F^ (Q^)/^ir (Q^) of 0(1) and 0(0?;,) contributions to the pion elastic form factor vs. Q^ for two reasonable values 
of the — here by assumption constant — effective threshold Seft. 



1. Since Scff((5^) is bounded from above, the 0(1) contribution i^ (Q^) to the elastic form factor F-^iCf') of the pion 

behaves like F)^ (Q^) oc l/Q'^ for Q^ — ?► oo. We would like to emphasize that the decrease of the soft contribution to 
the pion elastic form factor like 1/Q^ is a direct consequence of perturbation theory and quark-hadron duality. 

2. Consequently, for large Q^ the pion elastic form factor Ft^^Q"^) is dominated by the 0{as) contribution Fj^ (Q^)- 
Now, in the holographic models of [111, Il2[ merely the soft contribution is considered: it behaves like l/Q^ for 

large Q^. This is only possible if the effective threshold Seff(Q^) rises with Q^. However, this immediately leads to the 
violation of the factorization theorem for the 0{as) contribution, which is governed by the same effective threshold. 
Consequently, we would like to emphasize that the findings of [11, 12] would imply that the QCD factorization theorem 
is violated. We thus conclude that the predictions of [lllllJI for the pion form factor seem to us improbable.^ 



B. The three-point function (VAV) and the P ^- 7 transition form factor 

Let us now consider the amplitude of two-photon production from the vacuum induced by the axial- vector current of 
nearly massless quarks of one flavour, j^^ = 97^759, with £1^2 denoting the photon polarization vectors: 



{7{qih{q2)\j;ix = 0)|0) = r^„/,(p|gi, g2)e?e^, P^qi+q2. 



(2.14) 



The amplitude T^q/j is obtained from the vacuum expectation value of the T-product of two vector and one axial- 
vector currents and will be referred to as the (VAV) amplitude. Vector-current conservation yields the following 
relations: 



Tf,apip\qi,q2)qi = 0, T^c,p{p\qi,q2)q2 = 0- 

The general decomposition of the amplitude may be written as 

Tfial3(j}\qiiq2) = -Pli<^aPqiq2iFo + {qi<^iiafiq2 - qia<^liqil3q2)iFl + {q2'^p.Paqi ~ q2P<^liq2aqi)iF2- 

The form factor Fq contains the contribution of the pseudoscalar meson of our interest j3!^ such that 

Fp^y.y{ql,ql) = --— Fo(p2 = m],\ql,ql). 
JP 

We also consider the transition amplitude of the pseudoscalar current operator qj5q: 

{liqihiq2)\qi5q\0) = eo,i3g,g,e'ie^F5{qlqlp^). 



(2.15) 
(2.16) 

(2.17) 
(2.18) 



^ A way out would be to assume different effective tliresholds for the 0(1) and the 0{as) contributions to the pion form factor. This seems, 
however, a rather artificial construction. 



The two-photon amphtude of the divergence of the axial current takes the form 



{l{qihiq2Wm = 6„^,,,,e?e^(p^Fo - qfF, - q'.F^) 



Pi 



(2.19) 



The case of our interest is q\ = 0, then the form factor Fi does not contribute to the divergence. In perturbation theory, 
the form factors Fq, F2, and F^ may be written in terms of their spectral representations in p^ (with g^ = gl = — Q^): 



F.{p\q') 



ds 



■X{s,q^ 



(2.20) 



4m^ 



To one-loop order, the spectral densities read [201-122 

A2(s,g2) = - 
A5(s,g') = - 
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= ^yl - Am? /i 



(2.21) 



Obviously, the absorptive parts A,; obey the classical equation of motion for the divergence of the axial current 

s Ao(s, q^) - q^ A2(.s, q^) = 2m Asis, q^). (2.22) 

The form factors then satisfy 



00 
p'Fo{p\q^)-q^F2{p\q^)^2mF^{p\q^)-'^ J dsAois,q^). 



(2.23) 



4m^ 



The last integral is equal to — l/27r, independently of the values of m and q^, and represents the axial anomaly [9|: 



p^Foip^q^)~q^F2{p\q^) = 2mF,{p^,q^) + ^. 



(2.24) 



In the chiral limit m — and for q^ — 0, the form factor Fq develops a pole related to a massless pseudoscalar meson 
|34| . The residue of this pole is again the axial-anomaly 1/2tt^. Notice however that the pole in the one-loop expression 
for Fq disappears as soon as one of the photons is virtual 20, 34, 35(] . 

As is clear from (|2.23l) , the anomaly represents the integral of Aq , the spectral density of the form factor Fq . Adler 
and Bardeen [lOj tell us that the anomaly is non-renormalized by multiloop corrections. The easiest realization of this 
property would have been just the vanishing of multiloop contributions to the spectral density Ao(s, q"^). An argument 
in favour of this possibility comes from explicit two-loop calculations :24] which report the non-renormalizability of the 
full {VAV) vertex to the two-loop accuracy. However, if so, the full form factor Fq is given by its one-loop expression. 
Then, this expression should develop the pion pole, known to be present in the full amplitude for any value of q^. 
But, obviously, this pole does not emerge in the one- loop expression for Fq if g^ 7^ 0! 

This requires that multiloop corrections to the form factor Fq (and, respectively, to its absorptive part Ao(s, q'^)) do 
not vanish. Then one may ask oneself how it may happen that the anomaly nevertheless remains non-renormalized by 
multiloop corrections? The only possible answer we see [33| is that the non-renormalization of the anomaly is reached 
due to some conspiral property of multiloop contributions to Ao(s,q^) forcing their integral to vanish (in fact, quite 
similar to the one-loop result: although the spectral density explicitly depends on m and g^, its integral is an m- and 
q^-independent constant). 

As is obvious from (|2.16p , the contribution of the light pseudoscalar constitutes a part of the form factor Fq . The 
Borel sum rule for the corresponding Lorentz structure reads {Q^ — —q^ > 0) 

/ dsexp(— sr)Ao(s, Q^) = —fpFp^{Q^) exp(— mpr) -I- contributions of excited states. (2.25) 

Exploiting the concept of duality, the contribution of the excited states is assumed to be dual to the high-energy region 
of the diagrams of perturbation theory above an effective threshold Scff • After that, setting the Borel parameter t = 



we arrive at the LD sum rule for a pseudoscalar qq-raeson 

Sotf(Q') 



dsAo{s,Q') = -fpFp^iQ'). (2.26) 

The spectral density Aq to one-loop order is given by (I2.2ip : two- loop corrections were found to be absent [2J]. As 
discussed above, higher-loop corrections to Aq cannot vanish; so the l.h.s. of (|2.26p is known to 0{al) accuracy. 
In the chiral limit, the LD expression for the form factor for the one-flavour case is particularly simple (cf. [18ll43|): 

We would like to emphasize that the effective threshold in (|2.27p depends on Q^ . Apart from neglecting a1 and 
higher-order corrections to the spectral density Ag, no approximations have been done up to now: we have just 
considered the LD limit r = 0; for an appropriate choice of Scs{Q'^) the form factor may still be calculated exactly. 
Approximations come into the game when we consider a model for Seff((5^). 

Independently of the behaviour of s^siQ'^), the form factor at Q^ = is related to the axial anomaly: Fp^{Q) — 
l/(27r^/p). QCD factorization requires Scff(Q^) — >• 47r^/|, for large Q^ . So, the simplest model compatible with this 
requirement is obtained by setting Soff(Q^) = 47r^/p for all values of Q^ [l3]- This is the "conventional LD model" 
()1.6|) yielding for the neutral pion case the Brodsky-Lepage interpolating formula [1.21 ^ 

As already pointed out above, the conventional LD model for the effective continuum threshold (|1.6p . defined by 
the choice Scs{Q'^) — ScsiQ'^) = 47r^/^ for all Q^, or a slightly more sophisticated approach of Q are approximations 
which do not account for subtle details of the confinement dynamics. Consequently, it is important to understand 
the accuracy to be expected within this approach; in other words, it is important to obtain some reliable estimate of 
the expected deviations of the exact Seff(Q^) from its LD limit sld in the momentum region Q^ > 4-6 GeV^. To 
this end, in the next Section we take advantage of the fact that in quantum mechanics all the bound-state properties 
may be found exactly by solving the Schrodinger equation. On the other hand, also in quantum mechanics we may 
construct LD sum rules. 

3. EXACT VS. LD FORM FACTORS IN QUANTUM-MECHANICAL POTENTIAL MODELS 

The relevance and the expected accuracy of the LD model may be tested in those cases where the form factor 
Pp-yiQ^) is known, i.e., may be calculated by other theoretical approaches or measured experimentally. Then, the 
exact effective threshold may be reconstructed from (|2.27|) , in this way probing the accuracy of the LD model. 

To probe the accuracy of the LD model, we now consider a quantum-mechanical example: the corresponding 
form factors may be calculated using the solution of the Schrodinger equation and confronted with the results of the 
quantum-mechanical LD model, which is constructed precisely the same way as in QCD. For the elastic form factor, 
it is mandatory to consider a potential involving both the Coulomb and the confining parts; for the analysis of the 
transition form factor one may start with a purely confining potential. 

The basic object for quantum-mechanical LD sum rules is the analogue of the three-point correlator of field theory 

M 

r^^iE,E'^Q)^{r' ^0\j^Jiq)j^\r^O), Q ^ \q\. (3.1) 

Here, H is the Hamiltonian of the model; the current operator J{q) is determined by its kernel {r'\J{q)\r) = 
exp{iq ■ r) (5*^''''(r — r'). We do not take the spin of the current into account, therefore the basic quantum-mechanical 
Green function is the same for both types of form factors discussed above. 



^ In an alternative approach to the P-y form factor [36l . l37l | , the pseudoscalar meson is described by a set of distribution amplitudes of 
increasing twist which are treated as nonperturbative inputs. In our analysis, the deviation of the effective threshold Se£f(Q^) from its 
asymptotic value 'in'^fp corresponds to some extent to the contribution of higher- twist distribution amplitudes in the approach of [37j . 



A. Elastic form factor 

The elastic form factor of the ground state is given in terms of its wave function ^ by 



Fei(Q) = (*| J(q)|*) = /"d3rexp(iq.T)|*(r)|2= f d^k^{k)^'{k + q), Q = \q\. (3.2) 

Here, ^ is the ground state of the Hamiltonian 

fc2 a 
H=—-- + V,o^i{r), r^\r\. (3.3) 

Zm r 

Because of the presence of the Coulomb interaction in the potential, the asymptotic behaviour of the form factor at 
large values of Q is given by the factorization theorem |38|] 

i^ei(g)- > p^^^, i?3^|*(r = 0)p. (3.4) 

The quantum-mechanical LD sum rule for the form factor Fc\{Q) is rather similar to that in QCD: The double Borel 
transform {E -^ T, E' ^ T') of p.ip may be written in the form 



rN«(r, r, q)^ j dk' exp (^^^r^ J dk exp (-|^t^ ANR,(fc, k', Q) + r^,^_.(T, T', g). 



(3.5) 



where rpj|^(,^(T, T', Q) describes the contribution of the confining interaction and Ap(JJj(/c, fc', Q) is the double spectral 
density of Feynman diagrams of nonrelativistic perturbation theory, see Appendix A for details. 

Setting T' = T = leads to the LD sum rule, in which case rpj^j,^. vanishes [2^. The low-energy region of 
perturbative diagrams — below some effective continuum threshold kefi(Q) — is assumed to be dual to the ground-state 
contribution, which reads Rg Fe\{Q'^). Finally, we arrive at the following LD expression for the elastic form factor: 



'■g I ^"TpcrtV"-/ p 2 



R, = I dkp^^,ik) = ^ + am'^ + Oia'), (3.6) 



fcoff(Q) fectt(Q) 



^ei°(g) ^ ^ J dk J dk'A'^^,{k,k',Q). (3.7) 



^g 





The explicit result for Ap^^^{ki, ^2, Q) is given in Appendix A. 

The factorization formula (13.41) is reproduced by the LD sum rule p.7p if the momentum-dependent effective 
threshold behaves as 

kesiQ) > fcLD = (67r2 i?g)i/3. (3.8) 



B. Transition form factor 

The analogue of the tt"/ transition form factor in quantum mechanics is given by 

Ftr.nsiQ,E) = (Vl/|j(q)_l_|r = 0), (3.9) 

The case of one real and one virtual photon corresponds to £' = and Q ^ 0. At large Q, the transition form factor 
^trans(Q) = Ftrans{Q, E — 0) satisfics the factorizatiou theorem 

Fu^nsiQ) > ^^- 3.10 

Recall that the behaviour p.lOp does not require the Coulomb potential in the interaction and — in distinction to the 
factorization of the elastic form factor — emerges also for a purely confining interaction. 



The LD sum rule for the form factor i^trans(Q) is constructed on the basis of the same three-point function (|3.ip 
and has the form 



F. 



LD 



(Q) = 



1 



fcoft(Q) 



Rn 



I oo 

dkJdk'A^^,{k,k\Q). 



(3.11) 



Notice that the fc'-integration is not restricted to the low-energy region since we do not isolate the ground-state 
contribution in the initial state. The asymptotical behaviour p.lOp is correctly reproduced by Eq. (13.111) for 



kesiQ -> oo) = A;ld- 



(3.12) 



C. Quantum- mechanical LD model 



As is obvious from p.8|) and (|3.12p . the effective thresholds for the elastic and for the transition form factors have 
the same limit at large Q: 



kcffiQ ^-00)= kctfiQ ^ CX)) = fcLD- 



(3.13) 



The LD model emerges when one assumes that also for intermediate Q one may find a reasonable estimate for the 
form factors by setting 

kcs{Q) = ks{Q) = fcLD. (3.14) 

Similarly to QCD, the only property of the bound state which determines the form factor in the LD model is Rg. 

D. LD vs. exact effective threshold 

Let us now calculate the exact thresholds kcs{Q) and fcoff(Q) which reproduce the exact form factor by the LD 
expression; they are obtained by solving the LD sum rules p.7p and p. Ill) using the exact form factors on the left- 
hand sides of these equations. The deviation of the LD threshold fcLD from these exact thresholds measures the error 
induced by the approximation p. 141) and characterizes the accuracy of the LD model. 

For our numerical analysis we use parameter values relevant for hadron physics: m — 0.175 GeV for the reduced 
constituent light-quark mass and a = 0.3. We considered several confining potentials 

Konf(r)-a„(mr)", n- 2, 1,1/2, (3.15) 

and adapt the strengths (T„ in our confining interactions such that the Schrodinger equation yields for each potential 
the same value of the wave function at the origin, \I'(r = 0) = 0.078 GeV'^/^, which holds for CT2 — 0.71 GeV, ci = 0.96 
GeV, and ai/2 — lA GeV. The ground state then has a typical hadron size ^ 1 fm. 
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Fig. 3: Exact effective thresholds in quantum mechanics for the elastic (left) and the transition (right) form factors for different 
confining potentials. _Rg = \^{r = 0)| . 

Figure [3] presents the exact efi^ective thresholds. Independently of the details of the confining interaction, the 
accuracy of the LD approximation for the effective threshold and, respectively, the accuracy of the LD elastic form 
factor increases with Q in the region Q^ > 5 — 8 GeV^. For the transition form factor, the LD approximation works 
well starting with even smaller values of Q. 
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THE PION ELASTIC FORM FACTOR 



The 0(1) and 0{as) spectral densities are at our disposal, and the only missing ingredient for obtaining the form 
factor is s^s- We know the rigorous constraints on the effective threshold — at Q^ = from the Ward identity and at 
Q^ — >■ cx) from factorization — so it is easy to construct a model for Scff(Q^) by a smooth interpolation between these 
values. A simple parameterization with a single constant Qo fixed by fitting the data at Q^ = 1 GeV^ might read |39j 



Scff(Q^ 



4^V.^ 



1 



.(0)A 



1 + tanh 



Q^\ a.(0) 



Ql = 2.02 GeV . 



(4.1) 



According to Fig. [l] our interpolation perfectly describes the well-measured data in the range Q^ « 0.5-2.5 GeV^. 

Note that the effective continuum threshold Seff((5^) in Eq. (j4.ip approaches its limit sld already at Q^ « 4-5 
GeV^. For Q"^ > 4-5 GeV^, it practically coincides with the LD effective threshold of ^7^. Moreover, for Q^ > 5-6 
GeV^ the formula (|4.ip is pretty close to the model of Q. Obviously, the model labeled BLM in Fig. [I] provides a 
perfect description of the available Ft^{Q'^) data in the region Q^ = 1-2.5 GeV^. For Q^ > 3-4 GeV^, it reproduces 
well all the data, except for a point at Q^ = 10 GeV^, where it is off the present experimental value, which anyhow 
has a rather large error, by some two standard deviations.'^ 

Interestingly, in the region Q^ > 3-4 GeV^ the BLM model yields considerably lower predictions than the results 
of the different theoretical approaches presented in Refs. [Ill - [l4| . 

For a given result for the pion form factor, we define the equivalent effective threshold as the quantity which 
reproduces this result by Eq. (|1.3p . Figure U displays the equivalent effective thresholds recalculated from the data 
and from the theoretical predictions for the elastic form factor from Fig. [1] For the Ft^{Q'^) predictions of [lllll2l. [l4| . 
the corresponding equivalent effective thresholds Seff(Q^) recalculated from (|2.10p are depicted in Fig.|4l In all cases, 
they considerably exceed, for larger Q^, the LD limit sld dictated by factorization. Moreover, their deviation from 
Sld increases with Q^. 
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Fig. 4: Left: the "equivalent effective threshold" extracted from the data (red) vs. 
Right: equivalent thresholds for the theoretical predictions displayed in Fig. [T] 



the improved LD model (BLM) of [38 



The exact effective threshold extracted from the accurate data at low Q^ suggests that the LD limit may be reached 
already at relatively low values of Q^ w 4 — 8 GeV^. However, the results in the right plot imply that the accuracy of 
the LD model still does not increase — or even decreases — with Q^ even in the region Q^ ~ 20 GeV^ , in conflict with 
both our experience from quantum mechanics and the hint from the data at low Q^ . We look forward to the future 
accurate data expected from JLab in the range up to Q^ = 8 GeV^. 



^ It is virtually impossible to construct models compatible with all experimental results within Q'^ = 2.5-10 GeV'^, as revealed by closer 
inspection of Fig.[l] those approaches which hit the data at Q^ = 10 GeV^ overestimate the better-quality data points at Q^ ^ 2—4 GeV'^. 
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5. THE (7r°, ry,r?') -^ 77* TRANSITION FORM FACTORS 

A. {r],T]') -^ 77* 

Before discussing the 7r° case, let us consider the 77 and 77' decays. Here, one has to take properly into account the 
?7 — 77' mixing and the presence of two — strange and nonstrange — LD form factors. Following [5|, 140 . |41| we describe 
the flavor structure of 77 and 77' as^ 



1^) 



.uu + dd |_ , . , I A ,uu + dd „ 

-) cosip — |ss) smc/), \ri ) — \ — — ) sm0 + |ss) cos0, 0^39.3. 



V2 



V2 



(5.1) 



The 77 and 77' form factors then take the form 



Fr^jiQ^) = ^K^(g')cos0- i^,^(Q2)sin0, F^.^(Q2) 



—^Fn-yiQ'^) sin 4) + ^Fs-yiQ"^) cos < 



1 



(5.2) 

Here, Fnj{Q'^) and Fsj{Q'^) are the form factors describing the transition of the nonstrange and ss-components, 
respectively. The LD expressions for these quantities read 

F„,(Q2)^ 1 f dsa^;:lis,Q^), F,,(Q2) = 1 / d,a^l(s,Q2), 

Jn J Js J 



(5.3) 



4m? 



where Upert a-^id cr^J^^ denote cTpert with the corresponding quark propagating in the loop. Let us mention that for an 
isosinglet axial- vector current, a QCD axial anomaly contributes to the amplitude of interest [4^ and to the spectral 
densities of the form factors. This effect is of the order a^ and is not be expected to be important at large Q^. An 
overall agreement of the form factor from the LD model for rj and 77' mesons with the data speaks in favour of this 
expectation. 

In numerical calculations we set mu = rna = and rus = 100 MeV. The LD model involves two separate effective 
thresholds for the nonstrange and the strange components |4(| : 



„(") 



47rV, 



2 f2 



In ~ 1.07/,, 



.('') 



/. ~ 1.36/,. 



(5.4) 



According to the experience from quantum mechanics, the LD model may not perform well for small values of Q^, 
where the true effective threshold is smaller than the LD threshold; however, for larger Q^ the LD model in quantum 
mechanics gives accurate predictions for the form factors, as illustrated by Fig. |31 Figure [S] shows the corresponding 
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Fig. 5: LD predictions for 77 and 77' vs. experimental data from [44| (black) and [43| (red). 

predictions for 77 and 77' mesons. One observes an overall agreement between the LD model and the data, meeting the 
expectation from quantum mechanics. 



For comparison with the form factors obtained in a scheme based on the octet-singlet mixing, we refer to [42 
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B. 



77 



Surprisingly, for the pion transition form factor, Fig. [6l one observes a clear disagreement between the results from 
the LD model and the BaBar data [i6| . Moreover — in evident conflict with the rj and r]' results and the experience 
from quantum mechanics — the data implies that the violations of LD increase with Q^ even in the region Q^ w 40 
GeV^! The effective threshold extracted from the BaBar data is compatible with a linear growing function of Q^ 
with no sign of approaching the LD limit. 



Q^F„(Q^) [GeV\ 




LD 



10 15 20 25 30 35 40 



s^„(Q\Ge\/] 



Q^ [Ge/] 




5 10 15 20 25 30 35 40 



Q^ [Gel/] 



Fig. 6: The LD 717 form factor vs. data from [4J] (black) and 116;] (red) and the corresponding equivalent effective threshold. 



The P — ?► 7 transition form factors have been addressed recently in many publications (see, e.g., [36l l37l. l42l |46| - |54| | ) . 
However, no convincing explanation of the full picture of the P — >■ 7 transition form factors has been offered. In fact, 
it is hard to find a convincing answer to the question why nonstrange components in 77, 77', on the one hand, and in 
7r°, on the other hand, should behave so much differently? 



6. SUMMARY 

We presented the analysis of the pion elastic and the 7r°, ry, 77' transition form factors from the LD version of QCD 
sum rules. The main emphasis was laid on the attempt to probe the accuracy of this approximate method and the 
reliability of its predictions. Our main conclusions are as follows: 

• The elastic form factor: Our quantum-mechanical analysis suggests that the LD model should work increas- 
ingly well in the region Q^ > 4 — 8 GeV^, independently of the details of the confining interaction. For arbitrary 
confining interaction, the LD model gives very accurate results for Q^ > 20 — 30 GeV^. The accurate data 
on the pion form factor at small momentum transfers indicate that the LD limit for the effective threshold, 
s^fP = 47r^/^, may be reached already at relatively low values Q^ = 5 — 6 GeV^; thus, large deviations from the 
LD limit at Q^ = 20 — 50 GeV^ reported in some recent publications |ll| - [l4| appear to us rather unlikely. 

• The P — >■ 77* transition form factor: We conclude from the quantum-mechanical analysis that the LD 
model should work well in the region Q^ > a few GeV^. Indeed, for the -q — )■ 77* and 77' — >■ 77* form factors, 
the predictions from LD model in QCD work reasonably well. Surprisingly, for the tt — >■ 77* form factor the 
present BaBar data indicate an increasing violation of local duality, corresponding to a linearly rising effective 
threshold, even at Q^ as large as 40 GeV^. This puzzle has so far no compelling theoretical explanation. Our 
conclusion agrees with the findings of [33, |4y, |47| obtained from other theoretical approaches. 
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Appendix A: Perturbative expansion of Green functions in quantum mechanics 

We construct the perturbative expansions of both polarization operator and vertex function in quantum mechanics. 



1. Polarization operator 

The polarization operator 11(E) is defined by [2711 

U{E) = {r' = 0\G{E)\r^O), (A.l) 

where G{E) is the full Green function, i.e., G{E) = {H — E)^^ , defined by the model Hamiltonian under consideration 



H = Ho + Vir) 



Ho 



1f_ 
2to' 



(A.2) 



The expansion of the full Green function G{E) in powers of the interaction potential V has the well-known form 

G{E) = Go{E) - Go{E)VGo{E) + Go{E)VGo{E)VGo{E) + ■■■ , (A.3) 

with Gq{E) = {Hq — E)~^ . It generates the corresponding expansion of Ii-{E): 

n(£;) = Ho (£;) + Hi (£;) + ••• . (a.4) 

Explicitly, one finds 

1 f d^k 



no(^) 

^i{E) 



{27r)3 J |i-S' 



1 



d^fc d^k' 



;y((fc-fc') 



l\2\ 



{2^f J k^-E^-E 
We consider interaction potentials V{r) which consist of a Coulombic and a confining part: 



(A.5) 
(A.6) 

(A.7) 



Then the expansion (|A.4I) becomes a double expansion in powers of the Coulomb coupling a and the confining potential 
Fcont (see Fig. [T]). 




+ 



E E 



+ 




no(E) ni(E) 

Fig. 7: Expansion of the polarization operator in terms of Coulomb (wavy line) and confining (dashed line) interaction potentials. 

The contribution to n(_B) arising from the Coulombic potential is referred to as the perturbative contribution, Ilpcrt- 
The contributions involving the confining potential Vconf (including the mixed terms receiving contributions from both 
confining and Coulomb parts) are referred to as the power corrections, Hpowor- For instance, the first-order perturbative 
contribution reads 

d^fc d-^fc' Ana am f d'^k 



u[''He) 



1 



(27r)e 



fc2 



E^-E(k- k'f 87r2 

2m ^ ' 



fc2 



^ Ifel 



(A., 



2m 2iTi \ 2m 

The integral diverges but becomes convergent after applying the Borel transformation l/(a — E") — > cxp(— aT).^ The 
Borel-transformed polarization operator n(T) has the form [55J 



n(T) = npert(T) + npowcr(r), npe,.t(r) = i^^-^j l + V2^mTa + -rmr^Ta^ + 0{a'' 



(A.9) 



^ Note that the Borel transform of the Green function (H — E) ^ yields the quantum- mechanical time-evolution operator in imaginary time 
U{T) = exp{-HT). 
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In the LD limit, that is, for T — >■ 0, only npcrt(T') will be relevant. Nevertheless, as an illustration we provide also the 
result for the power corrections npowcr(r) for the case of a harmonic-oscillator confining potential T^onf (^) = muP'r^ j^: 



-L^powciA-' ) 1 



2-nT 



Z/2 



--uj^T^ (l + -V2^^^a] + —Lo^T^ 
4 V 12 / 480 



(A.IO) 



Let us point out that npower(T' = 0) vanishes, similar to QCD. The radiative corrections in npcrt(T') have a less singular 
behaviour compared to the free Green function, so the system behaves as quasi-free system. In QCD such a behaviour, 
frequently regarded as an indication of asymptotic freedom, occurs due to the running of the strong coupling as and its 
vanishing at small distances. Interestingly, in the nonrelativistic potential model this feature is built-in automatically. 
Now, according to the standard procedures of the method of sum rules, the dual correlator is obtained by applying a 
low-energy cut at some threshold fccff in the spectral representation for the perturbative contribution to the correlator; 



1 f ..,,. / k^ 



^duai{T,kci{) ^ —^ / dfcfc exp(-— T 



1 



mna (Trma)^ 



k 



3fc2 



+ 0(a3) 



n 



power 



(T). 



(A.ll) 



By construction, the dual correlator IlduaK^'j ^off ) is related to the ground-state contribution by 

nduai(T, fceff) = ng(r) = R^ exp{-E^T), R^ = |^g(r = o)p. 



(A.12) 



As we have shown in our previous studies of potential models, the effective continuum threshold defined according to 
(|A.12I) is a function of the Borel time parameter T. For T = 0, one finds 



n 



1 -, am -n 

dual (fccff , T = 0) = ^ fc^g + — fceff 



(A.13) 



2. Vertex function 



We now calculate the vertex function T{E, E', Q), defined by 

T{E, E',Q) = (r' = 0\G{E)J{q)G{E')\r = 0), Q=\q\ 



(A. 14) 



where J{q) is the operator which adds a momentum q to the interacting constituent. The expansions (|A.3p of the full 
Green functions G{E) and G{E') in powers of the interaction entail a corresponding expansion oiT{E, E' . Q), cf. Fig. [51 



T{E, E', Q) = Tq{E, E', Q) + Ti{E, E' , Q) 



(A.15) 





E E E' 



^ + 




E E' E' 



E E' 
ro(E,E',q) r^(E,E',q) 

Fig. 8: Nonrelativistic Feynman diagrams representing the lowest perturbative contributions to the vertex function F(_B, E' , Q). 



For the vertex functions Ti{E, £", Q), i = 0, 1, . . . , in (jA.lSp . their double spectral representations may be written as 

dz dz' . , . „, ^^^^g^ 



T,{E,E',Q)^ J 



— E ^ ~ E' 

2m -^ 2m ^ 



-A,{z,z',Q). 



The vertex functions Fi(_E', E' , Q = 0) and the polarization operators lli{E) satisfy the Ward identities 

n,(s)-n,(s') 



r,{E,E\Q^O) 



E-E' 



(A.17) 



15 



which are equivalent to the following relations between the corresponding spectral densities: 

lim Ai(z, z', Q) = S{z ~ z) pi{z). 



(A.18) 



We shall c onside r the double Borel transform ^ ^ T and £" ^ T': (a-Ey^ ~> exp(-aT), (a'-E')-'^ -^ cxp(-a'T') 
Equation (|A.18p leads to the following Ward identities for the Borel images: 



r,(T,T',Q = 0) = n,(T + T'). 



(A.19) 



a. One-loop contribution To to the vertex function 
The zero-order one-loop term has the form 



^oiE,E',Q) = ^J- 



d^k 



\ 2m / \ 2m J 

This may be written as the double spectral representation 

ToiE,E',Q)^f^' ^' 

J 9^ 

where 



2m E£-E 



-Ao{z,z',Q), 



Aoiz, z',Q) ^ -^ J d'kSiz - k^) S{z' - {k - q) 



(A.20) 



(A.21) 



--(^^^((^' — Q^)^-4.Q^<0). (A.22) 

Hereafter, we use the notations k = y/z and k' = \/z^. In terms of the variables k and k' , the function takes the form 

e{{z' -z- Q2)2 _ 42Q2 < 0) = e{\k ~Q\<k' <k + Q). (A.23) 



b. Two-loop contribution Fi to the vertex function 

We consider here only corrections related to the Coulomb potential, since power corrections induced by the confining 
interaction vanish in the LD limit. The two-loop 0(a) correction receives two contributions and has the form 



Ti{E,E',Q) 



1 



d^k d3fc' 



'ina 



(2^^)'-^ €l~E^-E'ik-ik'-q)f 



1 



1 



2m 2m 



(A.24) 



Having in mind the subsequent application of a double Borel transformation in E and E' , it is convenient to represent 

Fi as a sum of two terms, Fi — Ti'' +r[' , with 



t[^\e,e\q)^^ f 

1 V > ^^J ^^5 J ,j^ 
\2m 



d^k' 



I 



d?k 



E') (1^^ -e)J (k- (k' - q)f [fc2 _ (fc' ~ qf] 
d^k f d^k' 



87r5 7 ^k^ ^ E) (^^ - E') J {{k + q) - k'Y [k'^ - [k + q) 



r?iE,E',Q).^J 



am f d^fc d^k' 



ML -^ E — 

2m 2m 



E^- E' [k + q-k'f 



+ 



{k' -q)^-k^ {k + qf -k 



(A.25) 



The double Borel transformation in E' -^ T and E' — > T' is now easily performed. 



.(a) 



.(a) 



Let us start with Tl . One integration in Tl may be performed, leading to 

d^fc' f d^k 



t[''\e,e',q) 



167r2 



2m 



-E'){^-E)\y\ J (^il^-E){^-E')\k\_ 



(A.26) 
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The first term corresponds to the contribution of the "left" two- loop diagram in Fig.|8l i.e., with the potential before the 
interaction with the current J{q), while the second term is represented by the "right" two- loop diagram in Fig.|Sl The 
corresponding double spectral densities have a form very similar to Aq: 

/:^[i{k,k\Q) = ^^i^dfc-QI < fc' < fc + Q)0(O < k)e{{) < k'), A[i{k,k',Q)^A[^l{k',k,Q). (A.27) 

Explicit calculations yield the following double spectral densities of the two contributions to Tl ' related to the "left" 
and "right" two-loop diagrams in Fig. [51 



(6),, . , „. am 1 



A[">{k,k',Q) = 



327r6 Qk 



log' 



k' ~Q + k 



k' -Q-k 



log' 



k' + Q + k 



k' + Q-k 



Aflik,k',Q) = Aflik',k,Q). (A.28) 



At Q = 0, r^i"^ (r, T', g = 0) satisfies the Ward identity, F^"^ (T, T', Q = 0) = n["'> {T + T), whereas F^^^ (T, T', Q = 0) 
vanishes: F^ ^ (T, T', Q = 0) = 0. For large Q and T, T' 7^ 0, F^ ^ (T, T', Q) assumes a factorizable form (see Eq. (|X25)) ): 

Tf\T,T',Q) ^ l^no(r)no(T'). (A.29) 

At the same time, both Fo(T, T', Q) and F'^"-^ (T, T', Q) are exponentially suppressed for large Q and T, T' 7^ 0. Hence, 
F-j^ determines the large-Q behaviour of the vertex function. 

c. Dual correlator 

The dual correlator FduaKT", T', Q) is constructed in a standard way, by application of a low-energy cut to the double 
spectral representation of the perturbative contribution to (JA.16P : 

fccff(Q,T) fcrff(Q,T') 

Fduai(T,r',g)= / dk2kexp(~^T] f dk' 2k' expf-^T'^ A{z, z',Q) +Tpo^,,{T,T',Q). (A.30) 


By construction, the dual correlator corresponds to the ground-state contribution exp{~EgT) exp{—EgT') Rg Fg{Q). 
In the LD limit T = and T' — 0, FpowcrlT", T' , Q) vanishes and the ground-state form factor Fg{Q) is related to the 
low-energy part of the perturbative contribution considered above: 

Kii(Q) KiiiQ) 

f dk2k I dk'2k'A{k,k',Q)=Fg{Q)Rg, (A.31) 



with A(fc, fc', Q) = Ao(fc, k', Q) + A[]] (fc, fc', Q) + A^"^ (fc, fc', Q) + Afl (fc, k', Q) + A[%{k, fc', Q) + 0{a^). 

In order to provide the correct normalization Fg{Q = 0) = 1 of the elastic form factor Fg{Q), the effective thresholds 
should be related to each other according to kes{Q = 0) = kcS', then the form factor is correctly normalized due to the 
Ward identity (jA.18|) satisfied by the spectral densities. 
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